Anisotropy effects on Rashba and topological insulator spin polarized surface states: a unified 

phenomenological description 



Emmanouil Frantzeskakis' Q and Marco Grioni' 

'institut de Physique de la Matiere Condensee, Ecole Polytechnique Federate de Lausanne (EPFL), CH-1015 Lausanne, Switzerland 

(Dated: December?, 2011) 

Spin polarized two-dimensional electronic states have been observed in metallic surface alloys with giant 
Rashba splitting and at the surface of topological insulators. We study the surface band structure of these 
systems, in a unified manner, by exploiting recent results of k • p theory. The model suggests a different 
way to address the effect of anisotropy in Rashba systems. Changes in the surface band structure of various 
Rashba compounds can be captured by a single effective parameter which quantifies the competition between 
the Rashba effect and the hexagonal warping of the constant energy contours. The same model provides a 
unified phenomenological description of the surface states belonging to materials with topologically trivial and 
non-trivial band structures. 



I. INTRODUCTION 

The spin degeneracy of electronic states can be lifted at sur- 
faces or interfaces by the combined effect of the spin-orbit 
(SO) interaction and an out-of-plane electric potential gradi- 
ent 1 1 1. The latter appears as an effective magnetic field in the 
rest frame of the electrons giving rise to the Rashba-Bychkov 
(RB) effect. Spin polarized low dimensional states have been 
experimentally observed in the surface states of metals with 
a large atomic SO coupling |2-5|. Such 2D electron gases 
(2DEG) can be described by a nearly-free electron (NFE) 
model where the RB SO interaction enters as a k-linear cor- 
rection term in the effective Hamiltonian: 

Hnfe-SO = Kkxf^y - kycr^) + (1) 

2ms 

where b is the Rashba parameter, a scaled effective elec- 
tron mass (i.e. mj — m* /It?) and the Pauli spin matrices. 

Similar spin polarized 2D states have been observed in 
a new class of quantum matter: the 3D topological insula- 
tors (TI) ll6l-nTl. TI and RB systems are characterized by a 
similar momentum-dependent helical spin arrangement lfT2l . 
Not surprisingly, Hnfe so as introduced in Eq. ([T}, has been 
used as an approximation to describe the low-energy, long- 
wavelength, surface band structure of the 3D TI materials 
lfT3l[T4ll . In this case, the linear term is dominant and b = /ivp, 
where Vp is the Fermi velocity at the Dirac point. The spin po- 
larized surface states described by a TI effective Hamiltonian 
must reflect the topologically non-trivial band structure of the 
bulk. The RB surface states have a similar origin in the SO in- 
teraction, but correspond to trivial bulk topologies and always 
give rise to pairs of constant energy contours encircling the 
spin degenerate points |7 15 1. Despite these differences, Hrb 
(for RB systems) and Hji (for TI materials) are essentially 
identical. 

Greatly different values of ms, and hence of the quadratic 
term in Eq. ([T]i cannot be the only explanation for the topolog- 
ical differences reflected in the surface states of the two kinds 
of systems. For instance, the overall band structure of RB- 
split Pb quantum well states with very heavy effective masses 
fl6< is clearly different from that of 3D TI materials JSHD. 



Even more interestingly, the surface state of Bi2Se3 has been 
fitted with an effective mass of only O.lSme, a value lower 
than those encountered in many RB systems I.14J . 

Here we investigate which parameters in Hnfe-so determine 
the topological differences reflected in the surface states of 
RB and TI materials with the same structural symmetry. To 
this end, anisotropy effects are very critical and arise naturally 
when the isotropic RB Hamiltonian [Eq. ([T]l] is modified to 
describe real materials. The Rashba parameter in Eq. ([T]l was 
initially associated with the surface potential gradient 1 1 1, but 
large quantitative discrepancies in the size and anisotropy of 
the band splitting were observed because Hnfe-so disregards 
the important influence of the ion cores [2 |. One must there- 
fore take into account the atomic structure and the correspond- 
ing electronic charge distribution, either explicitly by using a 
site-spin representation lfT7HT9l . or implicitly by modifying 
the effective Hnfe-so- A single parameter cannot capture the 
anisotropic band structure and different approaches have been 
proposed. DFT calculations have stressed the importance of 
the wavefunction asymmetry near the ion cores 1201 |2TI . while 
an independent NFE approach enhanced the standard model 
with an ad hoc in-plane potential gradient |22 1. 

Recent data on Bi- and Pb-based surface alloys on metal- 
lic and semiconducting substrates have provided an experi- 
mental counterpart to the theoretical studies |23-28j . Those 
giant RB systems exhibit both a very large RB splitting and 
clear anisotropies. The anisotropy manifests itself in the shape 
of the constant-energy (CE) contours and the existence of an 
out-of-plane component of the spin polarization ll22l l23l l29l . 
Similar RB-type quantum well states with significant splitting 
have been also observed to coexist with TI surface states on 
several TI compounds where adatoms induce a strong band 
bending |'30'-'34 1 . In the field of TI, Fu has recently performed 
a k • p a perturbation analysis and proposed a higher order 
correction to Hti which can explain the anisotropy observed 
in the surface states of Bi2Te3 ["351. In this case, the extra 
parameter captures the deviation of their momentum distri- 
butions from perfect circles (i.e. the warping of the surface 
states). 

The purpose of this work is twofold. On one hand, we pro- 
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vide a phenomenological description which verifies that the 
results of k • p theory can be applied not only to TI but also to 
RB materials of the same symmetry (R3m). The k • p Hamil- 
tonian can be used for a common description of RB systems 
whose band structure follows Eq. ([T]l and of those where it 
deviates significantly from the predictions of the classical RB 
effect. Band structure changes in RB materials can be accu- 
rately captured by the continuous tuning of a dimensionless 
material-dependent parameter which quantifies the interplay 
of SO and anisotropy effects. On the other hand, using the 
same analysis, we investigate the possibility of using this pa- 
rameter to discriminate the spin polarized surface states of 
compounds with non-trivial bulk topologies from those ex- 
hibiting a RB-type dispersion. Our results have been used to 
model the surface band structure of real systems with remark- 
able accuracy around the center of the surface Brillouin zone 
(SBZ). 

II. MODEL AND PARAMETERS 

As developed in Ref. [35] the k • p Hamiltonian has to re- 
spect the symmetry class of the studied materials. Momentum 
and spin must be invariant under the necessary symmetry op- 
erations. Therefore, for the R3m group the system Hamilto- 
nian must take the following form up to third order terms: 

H = Hnfe-so + ^((fcx + ikyf + {K - ikyf)(Tz (2) 

The cubic term introduces warping effects on the otherwise 
circular in-plane momentum distributions. Deviations of the 
resulting dispersion from a single NFE parabola due to SO 
coupling and hexagonal warping can be quantified by cor- 
responding wavevector or energy scales. The characteristic 
wavevector scale for the isotropic RB effect, fci = b-nis deter- 
mines the fc-shift from the parent spin degenerate band, while 
the energy offset Ei = {b^ ■ TOs)/2 marks the energy posi- 
tion of the band's extremum. Namely, for m* > the split 
bands have a minimum at E = —\Ei\, while for m* < 
they have a maximum at E = \Ei\ |[3] l23l . Using dimen- 
sional analysis, one can define a characteristic wavevector and 
energy scale introduced by the hexagonal warping on a NFE 
parabola. k2 = l/(3c • rris) and E2 = l/(54c^ • m^) mark 
the /c -position and energy of the local extremum for the band 
dispersion normal to a mirror plane. For m* > the band has 
a local maximum at energy E — \E2\, while for ni* < the 
band has a local minimum at energy E — — | £'2 1 • A sketch of 
the effect of different parameters on a model band structure is 
given in Fig. |3jb). One should bear in mind that fci (fc2) and 
El {E2) are defined for a negligible value of c ih). When nei- 
ther 6 nor c is small, the momentum and energy values of the 
RB minimum and the local maximum are only approximately 
given by fci, k2, Ei and £'2. Therefore, in Fig. [Sj^b), k'l, fcj, 
£■* and E2 have been used to denote the corresponding mo- 
mentum and energy positions of the extrema. In the following 
analysis, we will use fci, fc2, Ei and E2 as the fundamental 



momentum and energy scales. The deviation of the "asterisk" 
values with respect to the corresponding fundamental param- 
eters increases with the intermixing of the Hamiltonian terms 
preceded by b and c. 

El and E2 are competing energy scales, since the former 
favors the RB minimum for m* > and the latter a local 
maximum in the dispersion of the outer spin polarized branch 
normal to the mirror planes [see Fig s. [3| b) and[3]^c)]. The 
characteristic wavevector scale a = \/b/c introduced by Fu is 
equal to ^/ki ■ 3fc2 and might be alternatively used to charac- 
terize the effects of RB coupling and hexagonal warping f35j . 
As noted by the author, that study discarded the effect of the 
effective mass because it does not contribute significantly to 
the shape of the constant-energy (CE) contours. Here, how- 
ever, we will not limit our conclusions to a single material. 
Our aim is to comment on the possibility of a generic clas- 
sification scheme by comparing spin polarized surface states 
with different mg, b and c. For this reason, one has to fol- 
low the modifications of the complete dispersion around F 
and the effective mass is one of the key parameters. In sec- 
tions III and IV we provide evidence that the ratio i? = fci / fc2 
(called effective anisotropy ratio hereafter) gives us the op- 
portunity to describe the splitting and warping effects using a 
single dimensionless constant. 2D electronic states with the 
same effective anisotropy ratio exhibit identical band topolo- 
gies irrespective of any differences in the values of TOs, b and 
c. R is also related to the ratio of the characteristic energies 

In order to solve the model of Eq. ([2]), we used a repre- 
sentation where the basis states are \k and \k \.), and spin 
projections refer to the z axis. If c = 0, the corresponding 
'spin-up' and 'spin-down' eigenstates refer to a quantization 
axis eg — ez X k/fc, which is always perpendicular to k, i.e. 
along the tangential direction. The corresponding states are 
in this case 100% in-plane polarized, with a purely tangential 
spin polarization P - opposite on the two branches - rotat- 
ing around F. For non-zero c the quantization axis deviates 
from the tangential direction resulting in non-negligible radial 
(Piad) and out-of-plane (Pz) spin polarization components. 
Experimental and theoretical studies have nevertheless proven 
that the in-plane rotation of the spins is preserved; i.e. despite 
the presence of Prad and Pz, the two eigenstates always refer 
to states with opposite signs of tangential spin polarization, 
while the total in-plane component is always tangent to the 
wai-ped CE contour Ii2"2l l23l |59l . 



III. QUALITATIVE RESULTS 

Fig. [T] summarizes our results for an increasing effective 
anisotropy ratio from (a) to (d). CE contours have been deter- 
mined from the eigenvalues of Eq. (|2]i. The results are plotted 
in units of fci, while the energy of the contours is equal to Ei. 
The image sequence clearly shows that there is a continuous 
transition from the isotropic RB effect (a) to the anisotropic 
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FIG. 1 : (color online) Model CE contours aiE = Ei for RB systems 
of R3rn symmetry characterized by different effective anisotropy ra- 
tios, (a) R = 0,(b)R = 0.225, (c) R = 0.292, (d) R = 0.45. Blue 
(black) and red (grey) colors denote opposite signs of tangential spin 
polarization. The mirror plane is along ky. The axes are scaled in 
units of fci . 



case (c). The outer contour is affected by the R3m symme- 
try of the surface and develops hexagonal warping due to the 
combined effect of the threefold axis and the time-reversal 
symmetry. The inner state also experiences the in-plane sym- 
metry of the surface as it is reflected at higher energies [Fig. 
|2|b)]. Figs. |2|a) and|2jb) evidence that the warping for the 
outer (inner) contour has a maximum positive (negative) value 
normal to the mirror plane (i.e. along fcx)^ while it vanishes 
along a mirror plane of the system (fcy). The shape of the CE 
contours follows the symmetry requirements imposed by Eq. 

The relative in-plane orientation of the inner and outer 
momentum distributions with respect to the mirror planes is 
in full agreement with the conclusions of a NFE model for an 
anisotropic RB effect on a 2DEG with R3m symmetry [22|. 

For even larger R values the outer constant energy contour 
is not closed around the T point [Fig. [ijd)]. In order to ac- 
cess its shape at larger wavevectors, one may renormalize the 
value of b using a recently proposed fc-quadratic correction 
[i.e. bk = b{l + ak^)] 1351 . Fig. [2|c) corresponds to the 
anisotropy ratio of Fig. [TJd) with a ~ 2. The outer state gives 
rise to six elliptical contours which shrink at higher energies 
while the inner spin split branch expands [Fig. |2|d)]. These 
contours are reminiscent of the Fermi surfaces of Bi(lll), 
Sb(lll) and the TI compound Bi , SK l4l l6l [TOl l36l [37l l40l . 
All of these systems are known to have a R3iti symmetry and 
exhibit two spin polarized branches. In each case, the outer 
branch does not form a closed contour around F but rather 
the so-called teardrop contours along the six equivalent mir- 
ror planes of the system. The teardrop states possess an un- 
conventional spin polarization whose in-plane component is 
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FIG. 2: (color online) [(a), (b)] Model CE contours at various en- 
ergy values for (a) the outer (E = 0.1 — 1.1-Ei) and (b) inner 
{E = 0.5 — 3.5-E2) branches of the system described in Fig. [T|c)- 
The contours expand at higher energies where their sixfold symmetry 
becomes more pronounced. The energy scales of (a) and (b) differ by 
about a factor of 10. [(c), (d)] The momentum distributions of Fig. 
[TJd) are extended to larger wavevectors by using a k-quadratic cor- 
rection for b. The results are plotted for E = E\ (c) and E = 2i?i 
(d) and resemble the measured Fermi surface of Bi(l 1 1) |4][36l and 
Sb(l 1 1) I37II . Blue (black) and red (grey) colors denote opposite val- 
ues of tangential spin polarization. The mirror plane is along fcy. The 
axes are scaled in units of fci . 



captured by the results presented in Fig. |2]([T0|. F is one of 
the time-reversal invariant momenta (TRIM) points. Accord- 
ing to theory an electron evolving by 360° in fc-space around 
a TRIM point acquires a geometric phase (Berry's phase) of tt 
BTI . Figs. [TJd),|2];c) and|2|d) suggest that an odd number of 
spin polarized contours encircle F, which might imply a total 
Berry's phase of tt as a consequence of the non-trivial topo- 
logical properties of the bulk. This is a necessary condition 
for the parent compounds of strong TI such as the Sb(lll) 
surface ifTOl 1411 . Although the present model cannot follow 
the surface state behavior at the other TRIM points of the sur- 
face Brillouin zone (SBZ) (i.e. M points), the band structure 
modification around F might be a hint that the value of R 
can discriminate between materials with topologically trivial 
and non-trivial band structures or consequently between pairs 
and an odd number of spin-polarized surface states around the 
TRIM points. Quantitative remarks and caveats about this de- 
scription will be given in the following. 

Fig. [3] compares the band dispersion for different effective 
anisotropy ratios and m* > 0. Along a miiTor plane direc- 
tion, the band structure remains unchanged because the warp- 
ing term of the Hamiltonian is odd and has no effect. Inter- 
esting changes can, however, be observed along the perpen- 
dicular direction: for increasing R, the dispersion becomes 
more anisotropic resembling the findings of Ref. 22 From 
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FIG. 3: (color online) [(a)-(d)] Model band dispersion for 2D states 
of R3rn symmetry characterized by different effective anisotropy ra- 
tios. The dispersion is plotted along and normal to a mirror plane. 

(a) R = 0, (b) R ^ 0.349, (c) R = 0.450, (d) R = 0.788. (e) 
The band structure of (d) is extended to larger wavevectors by in- 
cluding a fc-quadratic correction for b. The parameters introduced in 

(b) are explained in the text, (f) The energy value of the local band 
maximum as a function of R. The hatched area denotes the effective 
anisotropy ratios which are higher than R^^it (see text). Blue (black) 
and red (grey) colors in the band dispersion denote opposite signs of 
tangential spin polarization. The mirror plane is along ky. The axes 
are scaled in units of ki and Ei, while E — determines the energy 
position of the spin degeneracy point. 



(a) to (c), the outer branch continuously bends at lower en- 
ergy values and develops a local maximum [e.g. in Fig. [3] 
(c)] . When the values of the fundamental parameters E2 and 
El become comparable, the RB minimum becomes harder to 
define until it is finally not observed [Fig. |3jd)]. In order to 
access the band structure at larger wavevectors, we have plot- 
ted Fig. [3|d) by using the quadratic coiTection for b. The re- 
sults are presented in Fig. [3]^e). The dispersion is modified in 
a similar way along both high-symmetry directions by bend- 
ing the outer branch at even lower energy values. The effects 
are more visible at large values of k. As a result, the disper- 
sion resembles a warped cone as proposed by Fu ||33|. Un- 
der the same descriptive analysis, the shape of the electronic 
dispersion follows the suggestions of the electronic contours 
sequence presented in Fig. [T] 2D electronic states similar 
to those encountered in an anisotropic RB NFE gas may be 
considered as an intermediate case between the isotropic RB 
effect and surface states with a cone-like behavior accompa- 
nied by strong hexagonal waiping as observed in various TI 
compounds [10.,35.,39J . 

Fig. [Sj^f) presents a summary of the previous observations. 
It illustrates, for non-negligible b and c, the energy position 
of the local band maximum in a direction perpendicular to 
the miiTor plane, as a function of the effective anisotropy 
ratio. For increasing R, the energy position of the maxi- 
mum approaches the Rashba minimum at a decreasing rate 
(E ~ —El, due to the deviation of the fundamental parame- 
ters from the "asterisk" values as explained in section II). This 
is the lowest value it can have because no local maximum is 
observed for higher R values [Fig. |3jd)]. The energy value 
of the local maximum marks the transition from two closed 
concentric contours around T [e.g. Fig. [TJc)] to one closed 
contour and six additional elliptical pockets [e.g. Fig. [TJd)]. 
Therefore, the actual merit of Fig. |3jf) is that it determines the 
lower value of R (i?crit — 0.63) for which the spin polarized 
states always (i.e. at all energy values) exhibit an odd number 
of closed contours around F. 



IV. QUANTITATIVE RESULTS 

We will now use the results of the previous section in a 
quantitative analysis of published data on real systems. Table 
1 lists the effective anisotropy ratios, as well as the character- 
istic wavevector and energy scales for five different materials 
with R3in symmetry. Results of the fits are depicted in Fig. |4] 
and summarized in Table 1. The model of Eq. (j2]) has been 
solved using values of m* and b in good agreement with those 
in the corresponding Refs., while c is treated as the only free 
parameter 

The surface state of Au(ll 1) is included as the paradigm of 
an isotropic RB splitting. The consequences of the anisotropy 
start appealing for the spin split sp^ surface states of the 
BiCu2/Cu(lll) surface alloy which exhibits a non-zero R 
value. Nevertheless, anisotropy effects are not readily seen 
at the Fermi surface contours [Fig. |4|a) and Refs. [38]|42| be- 
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FIG. 4: (color online) Model band dispersion (upper panel) and Fermi surface contours (lower panel) with fit parameters which match published 
data (black dots) on real systems with R3fh symmetry. Black points have been digitally acquired from the Refs. mentioned in the last column 
of Table 1 . All parameter values of our model are summarized in Table 1 . (a) BiCu2/Cu( 111) spz surface state, (b) BiAg2/Ag( 111) sp^ surface 
state: The bottom panel refers to = —0.55 eV. Due to the lack of high-resolution experimental CE contours, black dots in the bottom panel 
refer to results from first-principles calculations [231. Before the fitting, the first-principles CE contours have been expanded to match the 
corresponding experimental data [i.e. the upper panel of Fig. [4jb) and Refs. 23 , 38 1. (c) Bi(l 1 1) surface state, (d) Cuo.o5Bi2Se3 surface state: 
The experimental dispersion along FK (top) has been obtained by Ref. 39 assuming a spherical symmetry up to Eb — —0.2 eV. The model 
results are compared to theoretical predictions (bottom), (e) (Bii_iSna)2Te3 surface state (5 = 0.67%): Corresponding fit parameters are in 
agreement with Ref. [35] In the bottom panel, the SBZ of the metallic surface alloys with a \/3 x \/3R30° surface reconstruction [i.e. (a), (b)] 
is rotated by 30° with respect to the other materials [i.e. (c)-(e)]. The dashed line denotes a mirror plane dictated by Eq. ([2j and it is along the 
FM direction of the SBZ. Blue (black) and red (grey) colors denote opposite signs of tangential spin polarization. 



TABLE I: Fit parameters of Eq. used to simulate (Fig. |4j the topography of different surface states with R3rh symmetry. Regular (bold) 
fonts denote surface states of RB (TI) materials. Parameter values follow the experimental dispersion published in the corresponding Refs. 
The energy value of the spin degeneracy point is set to zero. 



material 


m*(me) 


6(eVA) 


c(eVA^) 


fci(A-') 


fc2(A-') 


R 


-Bi(eV) 


S2(eV) 


-EF(eV) 


Refs 


Au(lll) 


0.27 


0.33 





0.012 


OO 





0.002 


OO 


0.475 




Bi/Cu(lll) 


-0.29 


0.85 


12 


-0.032 


-0.741 


0.043 


-0.014 


-2.439 


-0.215 


(Ml ED 


Bi/Ag(lll) 


-0.32 


2.95 


18 


-0.122 


-0.446 


0.274 


-0.181 


-0.800 


N/A 


l23l [2911381 


Bi(lll) 


0.50 


0.70 


58 


0.046 


0.088 


0.515 


0.016 


0.020 


0.067 


Elm 


Cuo.()5Bi2Se3 


0.26 


1.90 


140 


0.067 


0.068 


0.978 


0.063 


0.022 


0.45 


m 


(Bii_iSn«)2Te3 


1.37 


2.00 


230 


0.360 


0.008 


44.444 


0.360 


6x10"^ 


0.25 


011351 



cause the Fermi energy value is very small compared to £2- 
The corresponding spz states of the BiAg2/Ag(l 11) alloy are 
more anisotropic. This can be quantified by a substantial in- 
crease in the value of R. These surface states do not exhibit 
any Fermi surface contours but the anisotropy manifests itself 
in the CE maps taken at a binding energy Eg somewhat higher 
than the spin degenerate point [Fig. |4|b)]. 

The effective anisotropy ratio of the Bi(lll) surface states 
is even higher and Ep is now comparable to E2- The 
anisotropy effects are very strong on the Fermi suiface con- 
tours. They consist of a single spin polarized state centered at 
r and six hole pockets [4] |36| 1401 . which are reproduced by 
our results. As pointed out previously, such surface CE con- 



tours have been observed in materials with topologically non- 
trivial band structures ||6][T0l. According to the tight-binding 
model of Ref "T, the surface states of Bi(l 11) exhibit a com- 
plex dispersion at larger wavevectors. Therefore, the complete 
bandstructure can be revealed only by first-principles methods 
Nevertheless, our simple NFE-based results can provide 
valuable information about the dispersion in regions where the 
parent band can be approximated by a simple parabola, since 
they indicate that the outer branch forms contours which do 
not close around F. 

The 3D Tl compounds Bi2Se3 (Cu2;Bi2Se3 with x — 0.05) 
and Bi2Te3 (Sn-doped 0.27%) have been included as ex- 
amples of materials with R3m symmetry and high effective 
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anisotropy ratios. We have chosen these members of the TI 
family due to availability of high quality data which readily 
shows warping effects on the band dispersion and the CE con- 
tours |8, 39 1 . The fit parameters reproduce very well the ex- 
perimental results [Figs. |4|d) top and|4|e)] or the theoretical 
predictions [Fig. |4|d) bottom] for small fc-values. The outer 
spin polarized branch forms the lower part of the dispersion 
cone and there is a single spin polarized contour around the 
center of the SBZ, as expected for materials with a topologi- 
cally non-trivial band structure. The effective anisotropy ratio 
of Bi2Te3 is more than one order of magnitude higher than the 
one of Cu-intercalated Bi2Se3 as revealed by the correspond- 
ing parameters and reflected in the blossomlike shape of the 
Fermi surface contour. The hexagonal warping term c used 
for Bi2Te3 is in close agreement with the one proposed by Fu 
in his recent study ll35l . 

The prediction of an unmodified RB dispersion along the 
mirror plane direction of the 3D TI materials (i.e. not tak- 
ing into account a quadratic correction in the value of b) may 
seem at a first glance in striking contrast with the well-known 
cone-like dispersion of their surface states. Nevertheless, one 
should bear in mind that at energies slightly lower than the 
Dirac spin degeneracy point the surface states of TI materials 
merge into the bulk valence band continuum, thereby obscur- 
ing any clear experimental identification 181 [39l. Such a dif- 
ference in the dispersion between the mirror plane direction 
and its perpendicular has been clearly identified by ab initio 
calculations on various compounds of the TI family ||T3| . 

V. CONCLUSIONS 

Warped in-plane momentum distributions like those de- 
picted in Fig. |4]are a signature of a fc-dependent out-of-plane 
spin polarization component ll22l l23l l29l [35l . Accord- 
ing to previous theoretical studies and experimental results, 
vanishes along the three equivalent surface mirror planes 
ll22l|29l[35 l, one of which is indicated by the dashed lines in 
Fig. |4] while the other two are rotated by ±120°. The value 
of Pz follows a threefold symmetry determining the overall 
symmetry of the total spin polarization vector under the re- 
quirements of Eq. Q. Moreover, the orientation of the inner 
and outer CE contours is a fingerprint for the symmetry of the 
Pz vector ll22l|29]| . Therefore, the results presented in the bot- 
tom panel of Fig. |4]can qualitatively capture the fc-dependent 
out-of-plane rotation of the total spin polarization vector. 

As a tradeoff for its simplicity, the Hamiltonian proposed 
by Fu [i.e. Eq. ([2])] has a few shortcomings. Since it is de- 
rived by k • p perturbation theory, it is valid at wavevector 
values around F. In order to satisfy this need we did not ex- 
tend our results to other high-symmetry points of the SBZ. 
One should bear in mind that, similar to Ref. |35l we used 
scaled axes in Figs. [T|3] Scaling is necessary to maintain a 
universal character of the trends we have addressed. Never- 
theless, the validity of the corresponding results at points far 
from the SBZ center depends on the absolute fc-distance from 



zero. In this frame. Fig. |4]proves that the conclusions from the 
scaled-axes Figures can indeed be applied to real systems. An- 
other shortcoming of the model is related to the passage from 
anisotropic RB-type 2D electronic states to surface states ex- 
hibiting a dispersion similar to those of TI compounds. It is 
not straightforward to claim a topological change because our 
results are limited to the surface states and they do not con- 
sider the bulk valence and conduction bands of the described 
compounds. However, the energy position of the bulk bands 
continuum is crucial to determine whether there is a partner 
switching in the spin polarized surface state branches within 
the bulk gap and between two TRIM points |i4T| . The topo- 
logical differences between similar compounds may depend 
on bulk band characteristics such as a critical bulk band inver- 
sion which discriminates the trivial band structure of Bi(l 11) 
from the non-trivial band structure of Sb(l 1 1) 1 15 j. Although 
the described model cannot capture such bulk effects, the ef- 
fective anisotropy ratio is a measure of how fast the surface 
state dispersion deviates from the isotropic case. Therefore, 
it may determine the minimum size of the bulk gap around 
F which would allow a partner switching to occur. Although 
Eq. Q has already been used to describe the surface states of 
Bi2Te3 135^1, the derived dispersion along the direction which 
connects two TRIM points remains unchanged in the absence 
of the fc-quadratic correction in the value of b. As a matter 
of fact, the best indication for a possible topological change 
within the model is the predicted change in the total number 
of closed spin polarized contours around F. According to Fig. 
|3jf), there will always be an odd number of closed contours 
for R > 0.63. In any case, since a purely two-dimensional 
model is not enough to capture the complete band structure of 
a TI compound, our results on this topic rest as a descriptive 
analysis of the alleged topological change. Final predictions 
should be made after determining the surface state behavior at 
different TRIM points which cannot be captured by Eq. Q. 

Nevertheless, the value of the present work is not limited 
to a purely phenomenological analysis of the consequences of 
the aforementioned topological change reflected on the sur- 
face states. The description of the RB effect using the Fu 
Hamiltonian [Eq. Q], although overlooked in the past RB lit- 
erature, offers major advantages. As explicitly demonstrated 
by our qualitative and quantitative arguments, it presents a 
simple but instructive way to capture the effects of anisotropy 
on the surface band structure. Continuous tuning of the ef- 
fective anisotropy ratio permits a unified description of differ- 
ent RB systems with complex topologies without resorting to 
relativistic DFT calculations and deviating from the conven- 
tional RB description. The same Hamiltonian [i.e. Eq. (|2]l] 
can capture the band structure of (a) the classical RB effect 
as observed on Au(lll) and was originally accessed by the 
isotropic RB Hamiltonian [i.e. Eq. ([T]l] |2|, (b) the "weakly" 
anisotropic RB systems such as the metallic surface alloys 
which have been previously described by DFT calculations 
or by an NFE model with an ad hoc in-plane potential gra- 
dient 1221 l23l . and (c) "strongly" anisotropic cases which do 
not show a pair of closed contours around the spin degeneracy 
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points and have been previously described only by relativistic 
first-principles calculations |4|. Differences between the cor- 
responding RB surface states are accurately captured by the 
effective anisotropy ratio R, which can be determined either 
from an experimental CE contour or from the band dispersion 
perpendicular to the mirror plane direction. Alternatively, the 
simple R parameter can be determined by fitting relativistic 
first-principles calculations. 

Can we attribute a physical meaning to the effective 
anisotropy ratio? As already mentioned, R quantifies the 
warping effects on the surface state dispersion, which are al- 
ways accompanied by a non-zero out-of-plane spin polariza- 
tion component. The latter is intuitively linked to an in-plane 
electronic potential gradient. Premper et al. introduced ex- 
plicitly such an in-plane gradient in the standard RB model 
and found that it can both influence the splitting size and 
alter the circular shape of the surface state contours || 22l . 
There are however alternative views of the same phenomenon. 
Bihlmayer et al. proposed that it is the asymmetry of the sur- 
face state wavefunction which drives the size and anisotropy 
of the splitting [20). Recent first-principles calculations sug- 
gested that the asymmetric charge distribution can arise from 
the unquenched local orbital angular momentum in the pres- 
ence of SO coupling 143)1 . Because of their computational 
complexity, such approaches face practical limitations, and 
the effective anisotropy ratio R can be considered as an al- 
ternative easy-to-use effective parameter, which quantifies the 
warping effects induced to surface states of R3rii symmetry 
according to the predictions of the aforementioned physical 
models. Therefore, R is intimately related to an in-plane 
asymmetry of the electronic charge distribution, which in turn 
may correspond to an asymmetric surface state wavefunction 
or again to an anisotropic quenching of the local orbital angu- 
lar momentum. 



VI. SUMMARY 

We presented a phenomenological description of the simi- 
larities and differences between the surface states of RB and 
TI compounds by continuously tuning an anisotropy-related 
parameter We showed that the k • p Hamiltonian proposed by 
Fu ^35 1 can be used not only to describe the surface states of 
TI compounds but also as a tool to quantify anisotropy effects 
in RB systems. We propose a general classification scheme for 
2D spin polarized states of R3m symmetry. Materials which 
exhibit an isotropic RB splitting, and 3D topological insula- 
tors with a single Dirac cone on their surface can be consid- 
ered as the two opposite ends of our classification. Systems 
with an anisotropic RB splitting naturally fall into place as 
an intermediate category. This simple approach suggests a 
unified description of the RB surface states on different ma- 
terials which exhibit complex spin polarized contours without 
deviating from the simplicity of the conventional RB Hamilto- 
nian. We hope that it will motivate more fundamental efforts 
to compute the in-plane electronic charge distribution and the 



momentum and energy dependence of the spin polarization 
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